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Эволюционный и фрагментарный подходы 
к задаче о равномерной нагрузке 


В работе представлены результаты исследования задачи о равномерной нагрузке. Подобная задача 
возникает, в частности, при моделировании структурных элементов учебного процесса в высшем 
учебном заведении. Проблема заключается в распределении учебных дисциплин во времени таким 
образом, чтобы максимум нагрузки на студента был минимален. Предложены эволюционный подход 
для поиска оптимального решения и фрагментарный подход для построения допустимого решения 
задачи о равномерной нагрузке. Предложенные методы протестированы на наборе случайных задач. 


Введение 


Целью данной работы является исследование задачи о построении расписания, 
удовлетворяющего критерию минимума максимальной нагрузки. Исследуемая задача 
является подзадачей для общей задачи построения расписания занятий в высшем учеб- 
ном заведении. Различные постановки этой задачи рассматривались в русскоязычной 
литературе [1-3]. Известен в англоязычной литературе обзор [4]. Также из англоязычной 
литературы могут быть выделены работы [5-7]. В [8] предложена единая модель для 
задачи об университетском расписании предметов, основанном на учебной программе 
(Сиглсиит-Базе4 соигзе итеабБ те ргоМет). Общая задача является МР-полной [9], [10] 
в большинстве своих постановок. Исходя из этого для решения таких задач целесооб- 
разно рассматривать эвристические и метаэвристические алгоритмы. В [11] предлагается 
обзор различных метаэвристик для задач об университетском расписании (Ошуегзйу Я - 
тете ргоетз). Классом таких метаэвристик является класс эволюционных ал- 
горитмов, схема функционирования которых описывается в настоящей работе. В работе 
также предлагаются новые подходы к поиску приближенного решения для исследуемой 
задачи. 

Постановка задачи. Задача о равномерной нагрузке рассматривается далее как 
задача теории расписаний. Общая модель теории расписаний включает три конечных 
множества: множество машин М = (М, ,..„М м} ‚ множество работ 1 = {Л.Л} и постав- 


ленное во взаимное соответствие каждой работе множество операций О’ = {07,...,0/ } [12]. 
У 


В рассматриваемой задаче множество машин М будет иметь мощность 1. Будем по- 
лагать, что машина может выполнять одновременно более одной операции (единичным 
интервалом времени является одна неделя). Длина расписания К задана в условии, ею 
может быть длительность семестра или учебного года. Для каждой работы Л, Е] на 


множестве её операций О’ задан линейный порядок. Работы „Л, Е] предполагаются 
непрерывными, то есть для любых двух операций о0/ и 0/., разность времени начала 


5/., операции 0/,, и времени завершения /./ операции 07 равна нулю 5/.,-//=0. 
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Работа может быть начата в любой момент не позже А — [, +1, где /, — время выполне- 
ния работы „,. Следовательно, число возможных решений для данной задачи равно 
м=П&-ь+Ь: Под расписанием будем понимать множество моментов начала 
У 
выполнения работ $ = {5',..,5“1. В качестве критерия равномерности расписания 
будем рассматривать максимум объема учебной нагрузки в неделю на одного студента. 
Критерий может быть вычислен по формуле Г = тах (С, [,) ‚ где [, — текущий мак- 
р 
симум нагрузки, Г, — максимум нагрузки за период 1. [, вычисляется как Г, = [+ 
+р,‚› где Г, ‚, — нагрузка в период Г, полученная от / — 1 работ, р, — время выпол- 


нения работы „Л, в период 1. 


Далее предлагается ещё более сузить класс рассматриваемых задач с целью 
приближения к практической постановке. На протяжении семестра занятия по данной 
дисциплине могут повторяться с определённой периодичностью, например, одна 
пара в две недели, тогда при условии множество, что запланировано шесть занятий, 


длина (время выполнения) работы равна [, =р,`п,-К,, где р, - длина периода 


работы „Л,, п, — число повторений периода, А, — число крайних нулей, то есть таких, 


которые стоят в начале или в конце работы. Такая работа может быть представлена 
следующим вектором (2,0,2,0,2,0,2,0,2,0) ‚ если взять период (2,0). Работы, в которых 
последовательности операций с одинаковым временем выполнения и одинаковой 
позицией в последовательности периодично повторяются, в дальнейшем будут называ- 
ться периодическими работами. Каждая такая работа может быть задана указанием 
периода и числа его повторений. Поскольку при построении решения крайние нули 
не учитываются, то, например, работа с периодом (0,2) будет эквивалентна работе с 


периодом (2,0) при одинаковом числе повторений периодов в этих работах. 


1. Описание эволюционной модели 


Одним из современных направлений поиска оптимальных решений дискретных 
оптимизационных задач является подход основанных на использовании генетических 
(или более обще — эволюционных) алгоритмов. 

Этот подход приводит к понятию эволюционной модели, которая состоит из 
следующих элементов: 

— базовое множество решений Х, на котором ищется оптимальное решение; 

— оператор построения начальной популяции: оператор, который позволяет выделить 
на множестве Х его подмножество Хо Х; 

— оператор кроссовера К: ХхХ —>Х, позволяющий по двум перестановкам- 
родителям построить новую перестановку-потомок из множества Х; 

— оператор мутации М: Х > Х; 

— критерий отбора — алгоритм, позволяющий сравнивать по качеству решения в рамках 
заданной популяции; 

— оператор селекции, позволяющий выделять множество пар в популяции для выпол- 
нения кроссовера; 

— оператор отбора, позволяющий строить новые популяции из множеств родителей и 
потомков; 
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— условие остановки, которое определяет условие остановки эволюции. 
Эволюционный алгоритм поиска оптимального решения в эволюционной 
модели работает следующим образом. На каждом шаге с номером А формируется 
некоторое множество Л» < Х — текущая популяция. Для каждого элемента текущей 
популяции вычисляется значение критерия отбора. На начальном шаге Хо — это начальная 
популяция, построенная с помощью оператора построения начальной популяции. 

Далее с помощью оператора селекции выбирается множество пар текущей популя- 
ции для кроссовера. К каждой паре из выбранного множества пар применяются 
операторы кроссовера и мугации, которые позволяют получить множество потомков. 
Для каждого потомка также вычисляется значение критерия. С помощью оператора 
отбора на основе вычисленных значений критерия отбора из объединения текущей 
популяции и множества потомков формируется новая текущая популяция Аки, которая 
содержит элементы с наибольшими значениями критерия отбора. 

Работа алгоритма заканчивается при выполнении условия остановки. 


2. Представления решений для эволюционного алгоритма 


Использовались два представления решения задачи о равномерной нагрузке: 
векторами смещений и перестановками. 
Хромосома в виде вектора смещений (1....,й,...1,), 0<1<А-[, где К — длина 


расписания, [ — время выполнения работы „Л, может содержать такие 1, и 4, 


ТЕК, что Г, 


‚=&. Каждое 1 представляет разность между начальным моментом 


времени в расписании и моментом времени, в который начинается работа Л. 
Скрещивание хромосом осуществляется при помощи функции /(й,й), которая 
ставит в соответствие элементам вектора смещений родительских хромосом й,й, 
занимающим К-ю позицию в каждой из них, некоторое #›. Такой функцией могут быть 
тах(й,и), ши(й,) или (Е -й)-тт4 +1 , где т — генератор случайных чисел на 
интервале (0; ]. 
Оператор мутации заключается в случайной генерации номера позиции в 


векторе смещений и значения времени начала выполнения работы, соответствующей 
данной позиции. 

При представлении решения перестановкой элемент 7); перестановки Л пред- 
ставляет номер работы. Работы ставятся в расписание в порядке их вхождения в 
перестановку, при этом каждая работа ставится в расписание таким образом, чтобы 
минимизировать максимум нагрузки машины. Таким образом, на каждом шаге 
принимается локально-оптимальное решение. Оператор кроссовера работает таким 


образом, что получает на вход две родительские хромосомы-перестановки л' и л’, 
затем просматриваются элементы на первых ПОЗИЦИЯХ В ЭТИХ хромосомах Ив 
зависимости от значения 1-го элемента некоторого бинарного вектора (Ь.....,6,....,Ь, ) 


1-й ген хромосомы-потомка п? наследует значение гена из л', если Ь, =0, или из 


л?, если Ь; =1. Значение, унаследованное л?, удаляется из пил". 


Описанный кроссовер сохраняет свойство предшествования, то есть если ген 
со значением а предшествует гену со значением Ь в обеих родительских хромосомах, 
то данное отношение сохранится и в хромосоме-потомке. 
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Оператор мутации для перестановок генерирует две случайные позиции в 
перестановке и меняет их значения местами. 


3. Фрагментарная структура 


В соответствии с [13] фрагментарной структурой будем называть тройку (Х, У,К), 
где Х — конечное множество элементарных фрагментов, У — семейство подмножеств 
множества Х, К условие присоединения, т.е. правило, позволяющее ответить на 
вопрос, является ли объединение множеств из Х элементом множества У. 

Задача о равномерной нагрузке может быть сформулирована как задача с фраг- 
ментарной структурой. Фрагментами в этой задаче являются любые упорядоченные 
последовательности работ, соответственно в роли элементарных фрагментов выступают 
отдельные работы. Условие присоединения — это допустимое время начала присое- 
диняемой работы. Максимальный по включению фрагмент соответствует допустимому 
решению рассматриваемой оптимизационной задачи. Таким образом, допустимое реше- 
ние можно рассматривать как некоторую фиксированную перестановку работ л. 

В [13], [14] показано, как для отыскания допустимых решений подобных 
комбинаторных задач может быть использована фрагментарная структура задачи. 
Для построения допустимого решения применяется фрагментарный алгоритм Г. 
В общем случае, на вход фрагментарного алгоритма РЁ подаётся последовательность 
(1,1,,.... 1, ) номеров элементарных фрагментов, которая просматривается в поряд- 
ке нумерации, и фрагмент-решение 5, который в начальный момент пуст. На каждом 
шаге просмотра алгоритма ищется первый элемент в упорядоченной последователь- 
ности элементарных фрагментов, для которого выполняется условие присоединения. 
Это фрагмент добавляется к 5. Процедура повторяется до тех пор, пока не будет 
построен максимальный по включению фрагмент. 


4. Результаты тестирования 


Предложенные алгоритмы были реализованы на языке \У15иа! Ваз1с Юг АррйсаНоп$ 
и выполнены на процессоре Пе! Соте 2 Рио с тактовой частотой ядра 2 ГГц. Было 
сгенерировано 120 случайных задач. В этих задачах использовался период для работ 
вида (2,0). Для генерации использовался генератор задач, взятый из [15] с теми же 
начальными условиями. 

Результаты работы алгоритма сравнивались с работой метода случайного поиска. 
Параметры, которыми задавались индивидуальные задачи, были следующими: 

— А — длина расписания (10, 20, 50, 100); 
— п — число работ (10, 20, 50, 100); 
— 1,  — интервал, из которого случайно выбирается длина работы. 

Таким образом, получено 12 групп индивидуальных задач размерностью А на и 
по 10 индивидуальных задач в каждой группе. Результаты применения названных 
выше подходов к решению рассматриваемой задачи приводятся в табл. 1. 

Для каждой задачи каждый алгоритм запускался не менее 30 раз. Для полученных 
таким образом значений целевой функции было вычислено среднее. Значения целевой 
функции, которые приведены в табл. 1, представляют собой средние для каждой 
группы задач одинаковой размерности, вычисленные от средних для каждой индивиду- 
альной задачи. 

В наборе протестированных алгоритмов выделяются такие свойства: локальная 
оптимальность решения (ЛОР), наличие глобального поиска (в данном случае 
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эволюции). Эволюционный метод и случайный поиск и фрагментарный алгоритм с 
решениями, заданными на перестановках ЭА(пр), СП(пр), ФА, включают свойство 
ЛОР. В эволюционном алгоритме и случайном поиске с решениями, определенными 
векторами смещений ЭА(вс), СП(вс), ЛОР отсутствует. Для сравнения был отобран 
ФА с упорядочиванием по убыванию времени выполнения работ, поскольку именно 
с таким упорядочиванием ФА показал наилучшие результаты. 


Таблица 1 — Результаты применения различных подходов 


Алгоритм 
СП (вс) ЭА(вс) СП (пр) ЭА(пр) ФА 
Размерность 
10х10 8 8 8 7 7 
10х20 15 14 13 13 13 
10х50 40 36 33 32 32 
10х100 80 72 66 63 62 
20х20 16 15 13 13 13 
20х50 38 35 30 30 29 
20х100 76 70 62 61 58 
50х20 16 15 13 12 12 
50х50 38 35 30 30 28 
50х100 76 71 62 61 57 
100х50 38 35 31 30 29 
100х100 77 72 62 62 57 


Результаты эксперимента показали, что изменение свойств в алгоритме СП(вс) 
в следующем порядке ведет к улучшению качества решения: 

1. Эволюция. 

2. ЛОР. 

3. Эволюция и ЛОР. 

4. ЛОР и упорядочивание работ по убыванию их времени выполнения. 

Таким образом, можно сделать вывод о том, что эволюция и ЛОР улучшают по 
отдельности качество решения по сравнению со случайным поиском без ЛОР. При 
добавлении к ЛОР упорядочивания работ по убыванию их времени выполнения 
качество результатов повышается. 


Заключение 


В данной статье предложены два новых подхода к решению задачи о равно- 
мерной нагрузке. Получено представление решений в эволюционной и фрагментарной 
моделях для задачи теории расписаний. Выполнено тестирование предложенных 
алгоритмов на случайном наборе задач. 

Результаты тестирования, которые приведены в табл. 1, показывают перспек- 
тивность подхода на основе эволюционной и фрагментарной моделей для поиска 
приближенных решений задач теории расписаний. 
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О.С. Бондаренко, 1.В. Кон 

Еволющйний та фрагментарний шдходи до задач! про рвном не навантаження 

У робот! представлен! результати досллдження задач! про р1вномрне навантаження. Под1бна задача 
виникае, зокрема, при моделюванн! структурних елеменив навчального процесу у вищому навчальному 
заклад1. Проблема полягае у розподл! навчальних дисциплн у час! в такий спос1б, щоб максимальне 
навантаження на студента було м!нмальним. Запропоновано еволющйний шШдхд для пошуку 
оптимальних розв’язкв та фрагментарний шдхд до побудови допустимих розв’язюв. Запропонован! 
методи протестован! на набор! випадкових 1ндив1дуальних задач. 


0.5. ВопаагепКо, Г.И. Кот 

п Фе рарег Ше гези $ оРа зу ое ип огт 1оа@тз ргоМет аге ргезете4. А зппЙаг рго ет аг15ез, ш 
рагисчаг, \Теп то4деПте зётасвга! сотпропеп оЁа |еагитз ргосез$ ш а ишуегзИу. ТБе рго ет соп$15{5 
огаПосайпе оЁ соиг5ез ш Яте $0 фай з4еп’$ тахипа| юа@ те \аз шшииите4. Ап еуо[айопагу арргоась 
Юг ап орйта| зоНоп зеагсН ап4 Назтещагу арргоасВ Юг а Еаз1Ые зоайоп сопзёгасНоп аге ргорозе4. 
ТБе ргорозе4 те#од$ аге {е5%{е4 оп шзапсез оепегайе4 аё тапот. 


Статья поступила в редакцию 08.07.2009. 
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